skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Saito, T"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Nitrogen hydrides such as NH3 and N2H+ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N2H+ (1-0) across different environments of an external spiral galaxy, NGC 6946. We find a strong correlation (r > 0.98, p < 0.01) between the HCN (1-0) and N2H+ (1-0) intensities across the inner ∼8 kpc of the galaxy, at kiloparsec scales. This correlation is equally strong between the ratios N2H+ (1-0)/CO (1-0) and HCN (1-0)/CO (1-0), tracers of dense gas fractions (fdense). We measure an average intensity ratio of N2H+ (1-0)/HCN (1-0) = 0.15 ± 0.02 over our set of five IRAM-30m pointings. These trends are further supported by existing measurements for Galactic and extragalactic sources. This narrow distribution in the average ratio suggests that the observed systematic trends found in kiloparsec-scale extragalactic studies of fdense and the efficiency of dense gas (SFEdense) would not change if we employed N2H+ (1-0) as a more direct tracer of dense gas. At kiloparsec scales our results indicate that the HCN (1-0) emission can be used to predict the expected N2H+ (1-0) over those regions. Our results suggest that, even if HCN (1-0) and N2H+ (1-0) trace different density regimes within molecular clouds, subcloud differences average out at kiloparsec scales, yielding the two tracers proportional to each other. 
    more » « less
  2. We present new HCN and HCO+(J= 3–2) images of the nearby star-forming galaxies (SFGs) NGC 3351, NGC 3627, and NGC 4321. The observations, obtained with the Morita ALMA Compact Array, have a spatial resolution of ∼290–440 pc and resolve the innerRgal ≲ 0.6–1 kpc of the targets, as well as the southern bar end of NGC 3627. We complement this data set with publicly available images of lower excitation lines of HCN, HCO+, and CO and analyse the behaviour of a representative set of line ratios: HCN(3–2)/HCN(1–0), HCN(3–2)/HCO+(3–2), HCN(1–0)/CO(2–1), and HCN(3–2)/CO(2–1). Most of these ratios peak at the galaxy centres and decrease outwards. We compare the HCN and HCO+observations with a grid of one-phase, non-local thermodynamic equilibrium (non-LTE) radiative transfer models and find them compatible with models that predict subthermally excited and optically thick lines. We study the systematic variations of the line ratios across the targets as a function of the stellar surface density (Σstar), the intensity-weighted CO(2–1) (⟨ICO⟩), and the star formation rate surface density (ΣSFR). We find no apparent correlation with ΣSFR, but positive correlations with the other two parameters, which are stronger in the case of ⟨ICO⟩. The HCN/CO–⟨ICO⟩ relations show ≲0.3 dex galaxy-to-galaxy offsets, with HCN(3–2)/CO(2–1)–⟨ICO⟩ being ∼2 times steeper than HCN(1–0)/CO(2–1). In contrast, the HCN(3–2)/HCN(1–0)–⟨ICO⟩ relation exhibits a tighter alignment between galaxies. We conclude that the overall behaviour of the line ratios cannot be ascribed to variations in a single excitation parameter (e.g., density or temperature). 
    more » « less
  3. ABSTRACT OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source, which includes very high energy (VHE, $$E\gt $$ 100 GeV) $$\gamma$$-ray data taken by the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov telescopes) and H.E.S.S. (High Energy Stereoscopic System) imaging Cherenkov telescopes. The discovery of VHE $$\gamma$$-ray emission happened during a high state of $$\gamma$$-ray activity in July 2016, observed by many instruments from radio to VHE $$\gamma$$-rays. We identify four states of activity of the source, one of which includes VHE $$\gamma$$-ray emission. Variability in the VHE domain is found on daily time-scales. The intrinsic VHE spectrum can be described by a power law with index $$3.27\pm 0.44_{\rm stat}\pm 0.15_{\rm sys}$$ (MAGIC) and $$3.39\pm 0.58_{\rm stat}\pm 0.64_{\rm sys}$$ (H.E.S.S.) in the energy range of 55–300 and 120–500 GeV, respectively. The broadband emission cannot be successfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the data set well and a proton-synchrotron-dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the broad-line region to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be a flat spectrum radio quasar (FSRQ), in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL Lac and FSRQ objects. 
    more » « less
    Free, publicly-accessible full text available May 15, 2026
  4. Aims.Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei (AGN). We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. Methods.We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE;E > 100 GeV) gamma rays to radio with MAGIC, VERITAS,Fermi-LAT,RXTE,Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We characterized the variability by deriving fractional variabilities as well as power spectral densities (PSDs). In addition, we investigated images of the jet taken with VLBA and the correlation behavior among different energy bands. Results.Mrk 421 was in widely different states of activity throughout the campaign, ranging from a low-emission state to its highest VHE flux ever recorded. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions with indices around 1.5. We observe strong correlations between X-rays and VHE gamma rays at zero time lag with varying characteristics depending on the exact energy band. We also report a marginally significant (∼3σ) positive correlation between high-energy (HE;E > 100 MeV) gamma rays and the ultraviolet band. We detected marginally significant (∼3σ) correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study, hence indicating the exceptionality of this flaring event in comparison with the rest of the campaign. The 2010 violent activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these unprecedented radio features prevent us from firmly associating them to the specific flares recorded during the 2010 campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recentIXPEobservations. 
    more » « less
    Free, publicly-accessible full text available February 1, 2026
  5. The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very high-energy (VHE, 100 GeV <E< 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multi-wavelength data made modelling of the blazar’s spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: (i) a leptonic two-zone synchrotron-self-Compton scenario, and (ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E> 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, observed by MAGIC and theFermi-LAT telescopes, yielded a conservative 95% confidence upper limit ofz ≤ 0.244 for the redshift of this blazar. 
    more » « less
    Free, publicly-accessible full text available February 1, 2026
  6. Context.Blazars exhibit strong variability across the entire electromagnetic spectrum, including periods of high-flux states commonly known as flares. The physical mechanisms in blazar jets responsible for flares remain poorly understood to date. Aims.Our aim is to better understand the emission mechanisms during blazar flares using X-ray polarimetry and broadband observations from the archetypical TeV blazar Mrk 421, which can be studied with higher accuracy than other blazars that are dimmer and/or located farther away. Methods.We studied a flaring activity from December 2023 that was characterized from radio to very high-energy (VHE;E > 0.1 TeV) gamma rays with MAGIC,Fermi-LAT,Swift,XMM-Newton, and several optical and radio telescopes. These observations included, for the first time for a gamma-ray flare of a blazar, simultaneous X-ray polarization measurements with IXPE, in addition to optical and radio polarimetry data. We quantify the variability and correlations among the multi-band flux and polarization measurements, and describe the varying broadband emission within a theoretical scenario constrained by the polarization data. Results.We find substantial variability in both X-rays and VHE gamma rays throughout the campaign, with the highest VHE flux above 0.2 TeV occurring during the IXPE observing window, and exceeding twice the flux of the Crab Nebula. However, the VHE and X-ray spectra are on average softer, and the correlation between these two bands is weaker than those reported in the previous flares of Mrk 421. IXPE reveals an X-ray polarization degree significantly higher than that at radio and optical frequencies, similar to previous results for Mrk 421 and other high synchrotron peaked blazars. Differently to past observations, the X-ray polarization angle varies by ∼100° on timescales of days, and the polarization degree changes by more than a factor of 4. The highest X-ray polarization degree, analyzed in 12 h time intervals, reaches 26 ± 2%, around which an X-ray counter-clockwise hysteresis loop is measured withXMM-Newton. It suggests that the X-ray emission comes from particles close to the high-energy cutoff, hence possibly probing an extreme case of the Turbulent Extreme Multi-Zone model for which the chromatic trend in the polarization may be more pronounced than theoretically predicted. We model the broadband emission with a simplified stratified jet model throughout the flare. The polarization measurements imply an electron distribution in the X-ray emitting region with a very high minimum Lorentz factor ($$ \gamma\prime_{\mathrm{min}}\gtrsim10^4 $$), which is expected in electron-ion plasma, as well as a variation of the emitting region size of up to a factor of 3 during the flaring activity. We find no correlation between the fluxes and the evolution of the model parameters, which indicates a stochastic nature of the underlying physical mechanism that likely explains the lack of a tight X-ray/VHE correlation during this flaring activity. Such behavior would be expected in a highly turbulent electron-ion plasma crossing a shock front. 
    more » « less
    Free, publicly-accessible full text available March 1, 2026
  7. null (Ed.)
    ABSTRACT Both the CO(2–1) and CO(1–0) lines are used to trace the mass of molecular gas in galaxies. Translating the molecular gas mass estimates between studies using different lines requires a good understanding of the behaviour of the CO(2–1)-to-CO(1–0) ratio, R21. We compare new, high-quality CO(1–0) data from the IRAM 30-m EMIR MultiLine Probe of the ISM Regulating Galaxy Evolution survey to the latest available CO(2–1) maps from HERA CO-Line Extragalactic Survey, Physics at High Angular resolution in Nearby Galaxies-ALMA, and a new IRAM 30-m M51 Large Program. This allows us to measure R21 across the full star-forming disc of nine nearby, massive, star-forming spiral galaxies at 27 arcsec (∼1–2 kpc) resolution. We find an average R21 = 0.64 ± 0.09 when we take the luminosity-weighted mean of all individual galaxies. This result is consistent with the mean ratio for disc galaxies that we derive from single-pointing measurements in the literature, $$R_{\rm 21, lit}~=~0.59^{+0.18}_{-0.09}$$. The ratio shows weak radial variations compared to the point-to-point scatter in the data. In six out of nine targets, the central enhancement in R21 with respect to the galaxy-wide mean is of order of $${\sim}10{-}20{{\ \rm per\ cent}}$$. We estimate an azimuthal scatter of ∼20 per cent in R21 at fixed galactocentric radius but this measurement is limited by our comparatively coarse resolution of 1.5 kpc. We find mild correlations between R21 and carbon monoxide (CO) brightness temperature, infrared (IR) intensity, 70–160 µm ratio, and IR-to-CO ratio. All correlations indicate that R21 increases with gas surface density, star formation rate surface density, and the interstellar radiation field. 
    more » « less
  8. The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new 2 − 4″ (∼75 − 150 pc at 7.7 Mpc) observations at 2 and 3 mm covering the central 50″ (∼1.9 kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO + , HNC, CS, HC 3 N, N 2 H + , C 2 H, CH 3 OH, and H 2 CO. We complemented these with published 1 mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density Σ SFR on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 – the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO + , and HNC). We find that HCO + correlates best with Σ SFR , but the dense gas fraction ( f dense ) and star formation efficiency of the dense gas (SFE dense ) fits show different behaviours than expected from large-scale disc observations. The SBE has a higher Σ SFR , f dense , and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2−1)/CO(1−0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO + /HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies. 
    more » « less
  9. ABSTRACT The feedback from young stars (i.e. pre-supernova) is thought to play a crucial role in molecular cloud destruction. In this paper, we assess the feedback mechanisms acting within a sample of 5810 H ii regions identified from the PHANGS-MUSE survey of 19 nearby (<20 Mpc) star-forming, main-sequence spiral galaxies [log(M⋆/M⊙) = 9.4–11]. These optical spectroscopic maps are essential to constrain the physical properties of the H ii regions, which we use to investigate their internal pressure terms. We estimate the photoionized gas (Ptherm), direct radiation (Prad), and mechanical wind pressure (Pwind), which we compare to the confining pressure of their host environment (Pde). The H ii regions remain unresolved within our ∼50–100 pc resolution observations, so we place upper (Pmax) and lower (Pmin) limits on each of the pressures by using a minimum (i.e. clumpy structure) and maximum (i.e. smooth structure) size, respectively. We find that the Pmax measurements are broadly similar, and for Pmin the Ptherm is mildly dominant. We find that the majority of H ii regions are overpressured, Ptot/Pde = (Ptherm + Pwind + Prad)/Pde > 1, and expanding, yet there is a small sample of compact H ii regions with Ptot,max/Pde < 1 (∼1 per cent of the sample). These mostly reside in galaxy centres (Rgal < 1 kpc), or, specifically, environments of high gas surface density; log(Σgas/M⊙ pc−2) ∼ 2.5 (measured on kpc-scales). Lastly, we compare to a sample of literature measurements for Ptherm and Prad to investigate how dominant pressure term transitions over around 5 dex in spatial dynamic range and 10 dex in pressure. 
    more » « less